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New upper bounds on sphere packings I 

By Henry Cohn and Noam Elkies* 



Abstract 

We develop an analogue for sphere packing of the linear programming 
bounds for error-correcting codes, and use it to prove upper bounds for the 
density of sphere packings, which are the best bounds known at least for di- 
mensions 4 through 36. We conjecture that our approach can be used to solve 
the sphere packing problem in dimensions 8 and 24. 
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1. Introduction 



The sphere packing problem asks for the densest packing of spheres into 
Euclidean space. More precisely, what fraction of M. n can be covered by congru- 
ent balls that do not intersect except along their boundaries? This problem fits 
into a broad framework of packing problems, including error-correcting codes 
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and spherical codes. Linear programming bounds [D] are the most powerful 
known technique for producing upper bounds in such problems. In particular, 
[KL] uses this technique to prove the best bounds known for sphere packing 
density in high dimensions. However, [KL] does not study sphere packing di- 
rectly, but rather passes through the intermediate problem of spherical codes. 
In this paper, we develop linear programming bounds that apply directly to 
sphere packing, and study these bounds numerically to prove the best bounds 
known 1 for sphere packing in dimensions 4 through 36. In dimensions 8 and 24, 
our bounds are very close to the densities of the known packings: they are too 
high by factors of 1.000001 and 1.0007071 in dimensions 8 and 24, respectively. 
(The best bounds previously known were off by factors of 1.01216 and 1.27241.) 
We conjecture that our techniques can be used to prove sharp bounds in 8 and 
24 dimensions. 

The sphere packing problem in M. n is trivial for n = 1, and the answer has 
long been known for n = 2: the standard hexagonal packing is optimal. For 
n = 3, Hales [Ha] has proved that the obvious packing, known as the "face- 
centered cubic" packing (equivalently, the ^3 or D3 root lattice), is optimal, 
but his proof is long and difficult, and requires extensive computer calculation; 
as of December, 2002, it has not yet been published, but it is widely regarded 
as being likely to be correct. For n > 4 the problem remains unsolved. Upper 
and lower bounds on the density are known, but they differ by an exponential 
factor as n — ► 00. Each dimension seems to have its own peculiarities, and 
it does not seem likely that a single, simple construction will give the best 
packing in every dimension. 

We begin with some basic background on sphere packings; for more in- 
formation, see [CS]. Recall that a lattice in W 1 is a subgroup consisting of the 
integer linear combinations of a basis of M. n . One important way to create a 
sphere packing is to start with a lattice A C R n , and center the spheres at 
the points of A, with radius half the length of the shortest nonzero vectors in 
A. Such a packing is called a lattice packing. Not every sphere packing is a 
lattice packing, and in fact it is plausible that in all sufficiently large dimen- 
sions, there are packings denser than every lattice packing. However, many 
important examples in low dimensions are lattice packings. 

A more general notion than a lattice packing is a periodic packing. In 
periodic packings, the spheres are centered on the points in the union of finitely 
many translates of a lattice A. In other words, the packing is still periodic 
under translations by A, but spheres can occur anywhere in a fundamental 
parallelotope of A, not just at its corners (as in a lattice packing). 



1- W.-Y. Hsiang has recently announced a solution of the 8-dimensional sphere packing problem 
[Hs], but the details are not yet public. His methods are apparently quite different from ours. 
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The density A of a packing is defined to be the fraction of space covered 
by the balls in the packing. Density is not necessarily well-defined for patho- 
logical packings, but in those cases one can take a lim sup of the densities for 
increasingly large finite regions. One can prove that periodic packings come 
arbitrarily close to the greatest packing density, so when proving upper bounds 
it suffices to consider periodic packings. Clearly, density is well-defined for pe- 
riodic packings, so we will not need to worry about subtleties. See Appendix A 
for more details. 

For many purposes, it is more convenient to talk about the center den- 
sity S. It is the number of sphere-centers per unit volume, if unit spheres are 
used in the packing. Thus, 

W 2 „ 



(n/2)! 



since a unit sphere has volume 7r ra / 2 /(n/2)!. Of course, for odd n we interpret 
(n/2)! as T(n/2 + 1). 

In most dimensions, there are not even any plausible conjectures for the 
densest sphere packing. The only exceptions are low dimensions (up to perhaps 
8 or 10), and a handful of higher dimensions (such as 12, 16, and 24). The 
most striking examples are 8 and 24 dimensions. In those dimensions, the 
densest packings are undoubtedly the E% root lattice and the Leech lattice, 
respectively. The Eg lattice is easy to define. It consists of all points of K 8 
whose coordinates are either all integers or all halves of odd integers, and 
sum to an even integer. A more illuminating characterization is as follows: 
Eg is the unique lattice in R 8 of covolume 1 such that all vectors v in the 
lattice have even norm (v,v). Such a lattice is called an even unimodular 
lattice. Even unimodular lattices exist only in dimensions that are multiples 
of 8, and in R 8 there is only one, up to isometries of R 8 . The Leech lattice 
is harder to write down explicitly; see [CS] for a detailed treatment. It is the 
unique even unimodular lattice in R 24 with no vectors of length \/2. These two 
lattices have many remarkable properties and connections with other branches 
of mathematics, but so far these properties have not led to a proof that they are 
optimal sphere packings. We conjecture that our linear programming bounds 
can be used to prove optimality. 

If linear programming bounds can indeed be used to prove the optimality 
of these lattices, it would not come as a complete surprise, because other pack- 
ing problems in these dimensions can be solved similarly. The most famous 
example is the kissing problem: how many nonoverlapping unit balls can be 
arranged tangent to a given one? If we regard the points of tangency as a spher- 
ical code, the question becomes how many points can be placed on a sphere 
with no angles less than tt/3. Odlyzko and Sloane [OS] and Levenshtein [Lev] 
independently used linear programming bounds to solve the kissing problem in 
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Figure 1. Plot of log 2 S + n(24 — ra)/96 vs. dimension n. 

8 and 24 dimensions. (The solutions in dimensions 8 and 24 are obtained from 
the minimal nonzero vectors in the Eg and Leech lattices.) Because we know 
a priori that the answer must be an integer, any upper bound within less than 1 
of the truth would suffice. Remarkably, the linear programming bound gives 
the exact answer, with no need to take into account its integrality. By contrast, 
in most dimensions it gives a noninteger. The remarkable exactness seems to 
occur only in dimensions 1, 2, 8, and 24. We observe the same numerically in 
our case, but can prove it only for dimension 1. 

Figure 1 compares our results with the best packings known as of Decem- 
ber, 2002 (see Tables 1.1(a) and 1.1(b) of [CS, pp. xix, xx]), and the best upper 
bounds previously known in these dimensions (due to Rogers [Ro]). The graph 
was normalized for comparison with Figure 15 from [CS, p. 14]. 

2. Lattices, Fourier transforms, and Poisson summation 

Given a lattice Act", the dual lattice A* is defined by 

A* = {y | (x, )/)eZ for all x G A}; 

it is easily seen to be the lattice with basis given by the dual basis to any 
basis of A. The covolume |A| = vol(M ra /A) of a lattice A is the volume of any 
fundamental parallelotope. It satisfies |A||A*| = 1. Given any lattice A with 
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shortest nonzero vectors of length r, the density of the corresponding lattice 
packing is 

W 2 fr\ n 1 
JnJ^Xv lAf' 

and the center density is therefore (r/2) n /|A|. 

The Fourier transform of an L 1 function / : W 1 —> R will be defined by 

f(t) = [ f(x)e 27Ti ^ dx. 

PROPOSITION 2.1. Let a = n/2 - 1. If f :R n —>■ R is a radial function, 

then 

f(t) = 27r\t\- a / f(r)J a (27rr\t\)r n / 2 d ri 
Jo 

where "/( r )" denotes the common value of f on vectors of length r. 

For a proof, see Theorem 9.10.3 of [AAR]. Here J a denotes the Bessel 
function of order a. 

We will deal with functions / : IR n —> R to which the Poisson summation 
formula applies; i.e., for every lattice Act" and every vector v 6 R n , 

(2.i) E/^«)4E^ iM />). 

xeA II teA* 

with both sides converging absolutely. It is not hard to verify that the right- 
hand side of the Poisson summation formula is the Fourier series for the left- 
hand side (which is periodic under translations by elements of A) , but of course 
even when the sum on the left-hand side converges, some conditions are needed 
to make it equal its Fourier series. 

For our purposes, we need only the following sufficient condition: 

Definition 2.2. A function / : R n — ► R is admissible if there is a constant 
5 > such that \f(x)\ and \ f(x)\ are bounded above by a constant times 
(l + \x\)- n ~ 5 . 

Admissibility implies that / and / are continuous, and that both sides of 
(2.1) converge absolutely. These two conditions alone do not suffice for Poisson 
summation to hold, but admissibility does. For a proof for the integer lattice 
Z n , see Corollary 2.6 of Chapter VII of [SW]. The general case can be proved 
similarly, or derived by a linear change of variables. 

We could define admissibility more broadly, to include every function to 
which Poisson summation applies, but the restricted definition above appears 
to cover all the useful cases, and is more concrete. 



694 



HENRY COHN AND NOAM ELKIES 



3. Principal theorems 

Our principal result is the following theorem. It is similar in spirit to work 
of Siegel [S], but is capable of giving much better bounds. Gorbachev [Go] has 
independently discovered essentially the same result, with a slightly different 
proof. (He concentrates on deriving Levenshtein's bound using functions / 
for which / has fairly small support, but mentions that one could let the size 
of the support go to infinity.) 

Theorem 3.1. Suppose f : W l — > R is an admissible function, is not 
identically zero, and satisfies the following two conditions: 

(1) f{x) < for \x\ > 1, and 

(2) f(t) > for all t. 

Then the center density of n- dimensional sphere packings is bounded above by 

/(0) 
2"/(0)' 

Notice that because / is nonnegative and not identically zero, we have 
/(0) > 0. If /(0) = 0, then we treat /(0)//(0) as +oo, so the theorem is still 
true, although only vacuously. 

Proof. It is enough to prove this for periodic packings, since they come 
arbitrarily close to the greatest packing density (see Appendix A). In particu- 
lar, suppose we have a packing given by the translates of a lattice A by vectors 
vi,...,vn, whose differences are not in A. If we choose the scale so that the 
radius of the spheres in our packing is 1/2 (i.e., no two centers are closer than 
1 unit), then the center density is given by 

N 



2™|A| 

By the Poisson summation formula (2.1), 

xeA II teA* 

for all v € M. n . It follows that 

2 



E E/(* + u ;- w *) = tttE/(*) 

l<j,k<NxGA I ' teA* 



l<j<N 



Every term on the right is nonnegative, so the sum is bounded from below by 
the summand with t = 0, which equals iV 2 /(0)/|A|. On the left, the vector 
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x+Vj—Vk is the difference between two centers in the packing, so \x+Vj—Vk\ < 1 
if and only if x = and j = k. Whenever \x + Vj — Vk\ > 1, the corresponding 
term in the sum is nonpositive, so we get an upper bound of Nf(0) for the 
entire sum. Thus, 

N 2 f(0) 
Nf(0) > -^p, 

i.e., 

s< «°> 



2»/(0) 

as desired. □ 



This theorem was first proved by a more complicated argument, which is 
given in the companion paper [C]. 

The hypotheses and conclusion of Theorem 3.1 are invariant under rotat- 
ing the function /. Hence, we can assume without loss of generality that / has 
radial symmetry, since otherwise we can replace / with the average of its ro- 
tations. The Fourier transform maps radial functions to radial functions, and 
Proposition 2.1 gives us the corresponding one-dimensional integral transform. 

As an example of how to apply Theorem 3.1 in one dimension, consider 
the function (1 — |x|)x[-i,i](^)- It satisfies the hypotheses of Theorem 3.1 in 
dimension n = 1, because it is the convolution of X[-i/2,i/2] ( x ) with itself, and 
therefore its Fourier transform is 

/ sin7rA 2 

Thus, this function satisfies the hypotheses of Theorem 3.1. We get a bound 
of 1/2 for the center density in one dimension, which is a sharp bound. This 
example generalizes to higher dimensions by replacing X[-i/2,i/2] ( x ) with the 
characteristic function of a ball about the origin. However, the bound obtained 
is only the trivial bound (density can be no greater than 1), so we omit the 
details. In later sections we apply Theorem 3.1 to prove nontrivial bounds. 
It will be useful later to have the following alternative form of Theorem 3.1: 

Theorem 3.2. Suppose f : W a — ► M is an admissible function satisfying 
the following three conditions: 

(1) /(0) = /(0)>0, 

(2) f{x) < for \x\ > r, and 

(3) f{t) > for all t. 

Then the center density of sphere packings in W 1 is bounded above by (r/2) n . 
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Theorem 3.2 can be obtained either from rescaling the variables in The- 
orem 3.1 or from the following direct proof. For simplicity we deal only with 
the case of lattice packings, but as in the proof of Theorem 3.1 the argument 
extends to all periodic packings (and hence to all packings). 

Proof for lattice packings. For lattice packings, the density bound in the 
theorem statement simply amounts to the claim that every lattice of covolume 1 
contains a nonzero vector of length at most r. We will prove this first for lattices 
A of covolume 1 — e, and then let e — > 0+. For such lattices, 

E /(*) = i^b E /(*)> 

xeA teA* 

by Poisson summation. If all nonzero vectors in A had length greater than r, 
then all terms except /(0) on the left-hand side would be nonpositive. Because 
all terms on the right-hand side are nonnegative, we would have 

/(o) > £ f(x) = -L- £ f{t) > A. 

xeA ieA* 



However, 

/(0) _ /(0) 
l-e l-e 



> /(0), 



which is a contradiction. Thus, every lattice of covolume strictly less than 1 
must have a nonzero vector of length r or less, and it follows that the same 
holds for covolume 1. □ 

It seems natural to try to prove Theorem 3.2 by applying Poisson summa- 
tion directly to a lattice of covolume 1, but some sort of rescaling and limiting 
argument seems to be needed. We included the proof to illustrate how to do 
this. 

Logan [Lo] has studied the optimization problem from Theorem 3.2 in the 
one-dimensional case (for reasons unconnected to sphere packing), but we do 
not know of any previous study of the higher-dimensional cases. Unfortunately, 
these cases seem much more difficult than the one-dimensional case. 



4. Homogeneous spaces 

The space W 1 is a 2-point homogeneous space; i.e., its isometry group acts 
transitively on ordered pairs of points a given distance apart. By studying 
packing problems in homogeneous spaces, one can put Theorem 3.1 into a 
broader context, in which it can be seen to be analogous to previously known 
theorems about compact homogeneous spaces. 
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We start by reviewing the theory of compact homogeneous spaces. See 
Chapter 9 of [CS] for a more detailed treatment of this material. Suppose X 
is a compact 2-point homogeneous space. We assume that X is a connected 
Riemannian manifold, of positive dimension. We can write X as G/H, where 
(G,H) is a Gelfand pair of Lie groups. Then L 2 (X) is a Hilbert space direct 
sum of distinct irreducible representations of G, say 0^L o Xr For eacn h 
evaluation gives a map fj : X — > V*, because Vj turns out to consist of 
continuous functions. We define 

Kj(x,y) = {f j (x),f j (y)). 

This is a positive definite kernel: for every finite subset C C X, we have 

2 



E K j( x >y) 

x,y£C 



E /;(*) 



> 0. 



Because of G-invariance, Kj(x, y) depends only on the distance between x and 
y. This function of the distance is a zonal spherical function; we can define a 
way of measuring distance t(x,y) and an ordering of the Vj's so that Kj(x,y) 
is a polynomial Pj of degree j evaluated at t(x,y). In general, t maps X x X 
to [0, 1], and t(x, y) = 1 if and only if x = y (note that it is not a metric). For 
the unit sphere in M n , we take t(x,y) = (1 + (x,y})/2, and the polynomial Pj 
is the Jacobi polynomial Pj a,/3 \t), where a = f3 = (n — 3)/2. 

Now suppose C is a finite subset of X. We get inequalities on C from the 
fact that for each j, the sum J2x£C fj( x ) nas nonnegative norm. We can apply 
these inequalities as follows to get an upper bound for the size of C, in terms 
of the minimal distance between points of C: 

Theorem 4.1 (Delsarte [D]). Suppose 

m 
3=0 

with cij > for all j and f(t) < for < t < r. Ift(x, y) < r whenever x and 
y are distinct points of C, then 

\C\ < /(l)/a . 

Proof. Suppose C satisfies t(x,y) < r for all distinct x,y € C. Then 
consider 

x,yeC 

This sum is bounded above by |C|/(1) since t(x,y) < r unless x = y, and 
is bounded below by |C| 2 ao since / — ao is a positive definite kernel. Thus, 
\C\ < f(l)/a . □ 
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Theorem 4.1 is the analogue of Theorem 3.1 for compact homogeneous 
spaces. To see the analogy clearly, we need to study R n as a homogeneous 
space. 

We can write R n as G/H, where G is the group of isometries of K n and 
H = 0{n). Then we need to decompose L 2 (M ra ) in terms of irreducible repre- 
sentations of G. It is no longer a direct sum, but it can be written as a direct 
integral; specifically, L 2 (M n ) = f£°ir r dr, where 7i> is the irreducible represen- 
tation of G consisting of functions whose Fourier transforms are distributions 
with support on the sphere of radius r. 

We can find the zonal spherical functions as follows. The representation n r 
is generated by the functions x i— ► e 2m ^ x ' y ^ with \y\ = r, so ir* consists of 
functions on the sphere of radius r. The evaluation map from W 1 to ir* takes 
a point x £ M. n to the function y i— > e 2m ^ x ' y ^ on the sphere of radius r. Thus, 
the zonal spherical functions are given by 



(This of course depends only on \x\ — X2I, and can be evaluated explicitly in 
terms of Bessel functions using Proposition 2.1.) In other words, they are given 
by functions whose Fourier transforms are delta functions on spheres centered 
at the origin. See Section 4.15 of [DM] for a more detailed discussion of this 
point of view. 

Now the analogue of positive combinations of the zonal spherical func- 
tions Pj (t) from the compact case is radial functions with nonnegative Fourier 
transform, and we can see that Theorem 3.1 corresponds to 4.1. 



In one dimension, we have already seen how to use Theorem 3.1 to solve 
the (admittedly trivial) sphere packing problem. Based on numerical evidence 
and analogy with the kissing problem, we conjecture that it can also be used 
to get sharp bounds in dimensions 2, 8, and 24. For reasons to be explained 
shortly, it is more convenient to work with Theorem 3.2 instead of Theorem 3.1, 
so we shall do so; we can convert everything to the framework of Theorem 3.1 
by rescaling the variables. 

In each of dimensions 1, 2, 8, and 24, the densest known packing is a lattice 
packing, given by a lattice that is homothetic to its dual. This lattice is Z in 
dimension 1, the A2 root lattice (i.e., the hexagonal lattice) in dimension 2, the 
Eg root lattice in dimension 8, and the Leech lattice in dimension 24. See [CS] 
for information about these lattices. Each of these lattices except A2 actually 
equals its dual, but that is not true for A^. However, we can rescale A2 so that 
the rescaled lattice is isodual, i.e., isometric with its own dual (in this case, 
via a rotation). 




\y\=r 




5. Conditions for a sharp bound 
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Suppose A is any lattice of covolume 1, such as an isodual lattice, and / 
is a radial function giving a sharp bound on A via Theorem 3.2 (i.e., r is the 
length of the shortest nonzero vector of A). By Poisson summation, we have 

E /(*) = E /(*)• 

xeA x-eA* 

Given the inequalities on / and /, the only way this equation can hold is if / 
vanishes on A \ {0} and / vanishes on A* \ {0}. This puts strong constraints 
on / and /. When A is isodual, the vector lengths in A and A* are the same, 
so / and / must both vanish on A \ {0}. 

Of course, there are similar constraints on / for a sharp bound in Theo- 
rem 3.1 (as opposed to Theorem 3.2), but we prefer to work with this context, 
since the isodual normalizations are more pleasant, and are the standard nor- 
malizations for Eg and the Leech lattice. 

It is natural to try to guess / from our knowledge of its roots. For example, 
in one dimension we could try 

/(x)=a _^n(i-£)^(^)', 

which clearly satisfies f(x) < for \x\ > 1 and has the right zeros. In fact, one 
can compute its Fourier transform and check that / is nonnegative everywhere 
(it has support [—1, 1] and is positive in (—1, 1)), so it solves the sphere packing 
problem in dimension 1, in a different way from the function in the previous 
section. 

Unfortunately, it seems difficult to generalize this approach to higher di- 
mensions. One can generalize this function by replacing the sine function with 
a Bessel function (see Proposition 6.1), but that does not yield a sharp bound 
in dimensions greater than 1. Attempts to write down a product with zeros at 
the right places for a sharp bound lead to products that seem intractable. 

One important thing to note is that for a sharp bound above dimension 1, 
it is not possible for / to have compact support, as it does in the examples 
involving sine and Bessel functions. If it did, then / could not have sufficiently 
densely-spaced zeros. To be precise, if / is a radial function with support in 
the ball B(0, R) of radius R about the origin, then the common value f(r) on 
vectors of radius r satisfies 

f(r) = f f{t)e 2m ^ dt, 
Jb(o,r) 

where x is any vector with \x\ = 1. This defines an entire function of r, and 
for all complex r, 

| /(r) | < e 2*fl|r| f 

Jb(o,r) 
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so / is a function of exponential type, and Jensen's formula implies that / can 
have at most linearly spaced zeros (see Section 15.20 of [Ru]). However, the 
nonzero vectors in the Leech lattice have lengths y/2k for integers k > 1, and 
those in Eg have lengths y/2k for integers k > 0. The function / must vanish 
at those vector lengths, and these roots are too densely spaced for / to have 
compact support. Of course, / also cannot have compact support (because / 
vanishes on A* \ {0}). 

One might wonder whether the restriction to radial functions is misleading: 
perhaps a nonradial function could be constructed more naturally. We cannot 
rule out that possibility, but consider it unlikely. Even if / is not radial, a 
sharp bound implies that / and / must vanish on concentric spheres centered 
at the origin and passing through the nonzero points of A and A*, respectively. 
The simplest reason is that if / proves that A is optimal, then it proves the 
same for every rotation of A. Alternatively, after rotational symmetrization 
/ and / must vanish on these spheres, and the inequalities on their values 
then imply that they must have vanished before symmetrization (the average 
of nonnegative values vanishes if and only if the values all do). It would seem 
strange for / to vanish on these spheres without being radial, but of course we 
cannot rule it out. 



6. Stationary points 

We do not know how to use Theorem 3.1 to match the best density bound 
known in high dimensions, that of Kabatiansky and Levenshtein [KL]. How- 
ever, it provides a new proof of the second-best bound known, due to Leven- 
shtein [Lev]: 

An 

A < Jn/2 
~ (n/2)! 2 4™' 

where jt is the smallest positive zero of the Bessel function Jf. (For more 
information about the asymptotics of this bound and how it compares with 
other bounds, see page 19 of [CS], but note that equation (42) is missing the 
exponent in j™/ 2 -) We will show how to use a calculus of variations argument 
to find functions that prove that bound. This approach is analogous to that 
used by Levenshtein. Yudin [Y] has also given a proof of Levenshtein's bound 
that seems reminiscent of our general approach, but not identical. 

To construct a function / for use in Theorem 3.1, we begin by supposing 
that there is a function g such that f(x) = (1 — \x\ 2 )g(x) 2 , so that / auto- 
matically satisfies the inequality f{x) < for \x\ > 1. (We write g instead of 
g for convenience later.) Assume that g is radial, and has support in the ball 
of radius R about the origin; we discuss these assumptions later. Notice that 
nothing in our setup requires / to be nonnegative, so we must check for this 
property later. 
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We have / = g * (g + Lg), where * denotes convolution and 

1 11 B 2 
T = — — V — 

4vr 2 dt 2 ' 

3=1 J 

so that under the Fourier transform, L corresponds to multiplication by — \x\ 2 . 
We require /(0) = 1, i.e., / 8 „g = 1. We want to maximize /(0), subject to 
this constraint. Notice that 

/(0) = / g(0-t)(g + Lg) (t) dt = [ g(t) (g + Lg) (t) dt = f g(g + Lg), 

JR™ JR™ JR" 

because g is radial (and hence even). 

Perturb g to g+h, where h has integral zero (so that /(0) does not change) 
and has support in the ball of radius R. Then the first order change in /(0) is 

/ (gh + hg + gLh + hLg) = 2 ( h(g + Lg) 

JR n JR n 

(where the equality comes from integration by parts). 

In order to have this vanish whenever J Rn h = 0, the function g + Lg must 
be constant (within the support of g), so for some constant c we have 

g + Lg = cx R , 

where \r 1S the characteristic function of the ball of radius R about 0. Then 
it follows from Proposition 2.1 that 

(1 - \x\ 2 )g{x) = c(R/\x\) n / 2 J n/2 (27rR\x\). 

Now 2itR must be a zero of J ra / 2 for the right-hand side to vanish at 
| a; | = 1, and c is determined (given R) by 

l=j g=j g + Lg=( c Xr = cR n ^ 2 /{n/2)\ . 

JR n JR n JR n 

(The function g + Lg has the same integral as g, because integration by parts 
shows that Lg has integral 0.) 

In fact, we take 2ixR to be the first positive root j n / 2 of J n /2, in order 
to make / nonnegative everywhere. To check that it is indeed nonnegative 
everywhere then, we use the equation 



f = g*(g + Lg)=g*(cx 



R) 



If g is nonnegative everywhere, then so is /. We can now determine g explicitly: 
We know that g + Lg = c Xr . The differential operator l + L takes a radial 
function u(\t\) to 

u"(\t\) n-l ... .. .. .. 



702 



HENRY COHN AND NOAM ELKIES 



It follows from this and the differential equation satisfied by the Bessel func- 
tions that 

J a {2ir\t\)/\t\ a 

is in the kernel of 1 + L, if a = n/2 — 1. From that, one can check that 
,s ( Ra Ja(2ir\t\) A 

One can check that the minimum of J a (27r\t\)/\t\ a occurs at 

1*1 = J„/ 2 /(27r) = R, 

so g is nonnegative everywhere, as desired, since J a {2itR) = J a {j n /2) < (the 
roots of the functions J a and J Q +i = J n / 2 are interlaced). 

Thus, we have constructed a function / satisfying the hypotheses of The- 
orem 3.1. It has /(0) = 1 and 

(n/2)! 



/(0) = f g(g + Lg)=c f g = c 



7T n / 2 R n ' 



Because R = J ra /2/(2vr), this function proves Levenshtein's sphere packing 
bound: 

Proposition 6.1. The function 

,, x ^n/2(jn/2kl) 2 

/W (l-|x| 2 )|x|™ 
satisfies the hypotheses of Theorem 3.1, and leads to the upper bound 

An 

Jn/2 



(n/2)! 2 4 n 

for the densities of n- dimensional sphere packings. 

This function does not optimize the bound in Theorem 3.1, but it does 
optimize it within the class of functions whose Fourier transforms have sup- 
port in the ball of radius jVi/2/ 71 " about the origin. This was first proved by 
Gorbachev [Go]. For another proof, see [C]. 

It is not a coincidence that this proves exactly the same bound as in 
Levenshtein's paper [Lev]. Levenshtein studies spherical codes with minimal 
angular separation 9, and derives his sphere packing bound from letting 6^0. 
Under that limit, the functions he uses in the linear programming bounds for 
spherical codes become our Bessel function example. This can be proved using 
the limit 

(6.1) lim j-p^ ( cos -) = (z/2)- a J a (z), 

which is 10.8 (41) in [EMOT]. 
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The functions we have obtained are not optimal in any dimension above 1. 
There are two reasons for this. First, we restricted our attention to functions 
such that / has compact support, and as we have seen in Section 5, that cannot 
be true if we are to get sharp bounds. Second, and more importantly, we 
implicitly considered only functions such that / is positive within its support. 
The problem is that if / vanishes somewhere, then the perturbations h must 
be chosen so as not to push f below zero. We made no attempt to do so, and 
therefore could not find any stationary points for which / has zeros within its 
support. We have seen in Section 5 that zeros are essential for sharp bounds. 

Unfortunately, it seems difficult to adapt the stationary point argument 
to deal with these difficulties. One approach is to assume that / has zeros at 
certain locations, and look at only the perturbations h that, up to first order, 
do not push / below zero. Although we can set up such problems, we have not 
been able to solve them. 

7. Numerical results 

It is possible to get numerical results by using linear programming to 
find functions for use in Theorem 3.1, as was done for the kissing problem by 
Odlyzko and Sloane in [OS]. The idea is to fix one of /(0) and /(0), and view 
extremizing the other as an infinite-dimensional linear programming problem. 
One can try to approximate it with a finite-dimensional problem, and solve it 
on a computer. Although we obtained some numerical results this way, it was 
cumbersome and generally ineffective. Instead, we use the following approach. 

First, consider trying to use our techniques to bound the density of an 
isodual lattice. There is no reason for optimal sphere packings to be isodual 
lattices, and for example in three dimensions they are known not to be, but it 
is convenient to use this case as a stepping stone. 

Proposition 7.1. Suppose g : R n — > R is a radial, admissible function, 
is not identically zero, and satisfies the following three properties: 

(1) 9(0) = 0, 

(2) g(x) > for \x\ > r, and 

(3) g = -g. 

Then every isodual lattice in dimension n must contain a nonzero vector of 
length at most r. 

Proof of special case. For simplicity, we deal only with the case in which 
g(x) > for \x\ 3> 0. Let A = A* be an isodual lattice. By Poisson summation, 

9( x ) = Yl 9(x) = ~J2 9{x), 

xGA rreA* :reA 
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SO 

E 9{x) = 0. 

x<=A 

In order for the sum not to be positive, the lattice A must contain some nonzero 
vector of length at most r. □ 

We can actually remove the hypothesis that g(x) > for \x\ 3> from the 
proof of Proposition 7.1, by using a scaling trick, as in the proof of Theorem 3.2. 
However, we omit the details, because the hypothesis holds in all our numerical 
examples. 

Notice that given any function / that proves a bound in Theorem 3.2, 
we can produce a g that proves the same bound in Proposition 7.1, by taking 
g = f — f. Thus, the bound we get for isodual lattices is at least as good 
as that for general sphere packings. In principle, it could be better, but in 
practice we find that it is not (as we explain below). 

We can find excellent functions for use in Proposition 7.1 as follows. Let 
L^(x) be the Laguerre polynomials orthogonal with respect to the measure 
e~ x x a dx on [0, oo). Set a = n/2 — 1, and define 

9k {x) = Ll(2^\x\ 2 )e- n \ x \ 2 

for k > 0; we suppress the dependence on n in our notation. These func- 
tions form a basis for the radial eigenfunctions of the Fourier transform, with 
eigenvalues (— l) h (see Section 4.23 and equation (4.20.3) of [Leb]). 

To find a function g for use in Proposition 7.1, we consider a linear com- 
bination of <7i , </3 , . . . , <?4 m +3 j an d require it to have a root at and m double 
roots at zi, . . . , z m . (Counting degrees of freedom suggests that there should 
be a unique such function, up to scaling.) We then choose the locations of 
Zi,...,z m to minimize the value r of the last sign change of g. To make this 
choice, we do a computer search. Specifically, we make an initial guess for 
the locations of z\ , . . . , z m , and then see whether we can perturb them to de- 
crease r. We repeat the perturbations until reaching a local optimum. Strictly 
speaking, we cannot prove that it ever converges, or comes anywhere near the 
global optimum. However, it works well in practice. At any rate this cannot 
affect the validity of our bounds, only their optimality given m. As we increase 
m, this method should give better and better bounds, which should converge 
to the best bounds obtainable using Proposition 7.1. 

This method gives good functions for use in Proposition 7.1, but naturally 
bounds on all sphere packings would be better than bounds only on isodual 
lattices. We can turn these functions into functions satisfying the hypotheses 
of Theorem 3.2, without changing r, as follows. 

Let ft be a linear combination of go,g2, ■ ■ ■ ,g4m+2 with double zeros at 
zi,...,z m , such that g + h has a double zero at r. Then in the examples we 
have computed, g + h has constant sign, which we can take to be positive. We 
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end up with a function / = —g + h that is nonpositive outside radius r, and 
whose Fourier transform / = g + h is nonnegative everywhere; furthermore, 
/(0) = /(0). Thus, we can apply Theorem 3.2 to get the same bound for 
general sphere packings that g proves for isodual lattices. Notice that it is not 
clear a priori that / must satisfy all the hypotheses of Theorem 3.2, but it 
happens in all of our numerical examples. 

Figure 1 (from §1) and Table 3 (from Appendix C) illustrate the bounds 
this method produces, using m = 6 as the number of forced double zeros 
(except in dimension 7 and lower, where m = 5 suffices for the accuracy we 
desire). We also list for comparison Rogers' bound [Ro], which was previously 
the best bound known in dimensions 4 through 36. The choices of forced double 
roots are described in Table 4 (from Appendix C). 

These bounds were calculated using a computer. However, the mathe- 
matics behind the calculations is rigorous. In particular, we use exact rational 
arithmetic, and apply Sturm's theorem to count real roots and make sure we 
do not miss any sign changes. More precisely, we take the quantities 2-Kzf 
to be rational numbers. That is convenient, because the functions gk{x) are 
polynomials in 27r|x| 2 with rational coefficients (times Gaussians, which have 
no sign changes). The one subtlety is that in constructing h, we want g + h to 
have a double root at r, and 2irr 2 is generally not rational. Instead of doing 
this computation using floating point arithmetic, we replace 2irr 2 by a nearby 
rational number and then determine h exactly, using that value of r. To do so, 
we must increase r slightly, but of course we can make the increase as small as 
we wish. 

In dimensions 8 and 24, we carried out the calculations for m = 11. The 
resulting upper bounds are within factors of 1.000001 and 1.0007071 of equality, 
respectively. More precisely, in dimension 8 we take 

2vrr 2 = 12.56637375, 

and in dimension 24 we take 

27rr 2 = 25.1342216. 

The forced double roots we used to achieve these bounds are given in Tables 1 
and 2. 





1 


2 


3 


4 


5 


6 


2irzf 


25.1328 


37.6995 


50.2678 


62.8463 


75.4579 


88.2463 





7 8 


9 


10 


11 


2nzf 


101.3618 115.6443 


131.0298 


150.0861 


174.2876 



Table 1. Forced double roots for m = 11, n = 8. 
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1 


2 


3 


4 


5 


6 


2^ 2 


37.705 


50.285 


62.893 


75.578 


88.454 


101.737 





7 8 


9 


10 


11 


2nzf 


115.776 131.035 


148.162 


168.215 


193.766 



Table 2. Forced double roots for m = 11, n = 24. 

Our numerical results lead us to the following conjectures. 

Conjecture 7.2. The smallest possible value of r in Proposition 7.1 
equals that in Theorem 3.2, and for each optimal g from Proposition 7.1, there 
exists an optimal f from Theorem 3.2 such that g = f — /. 

Conjecture 7.3. There exist functions that satisfy the hypotheses of 
Theorem 3.2 and solve the sphere packing problem in dimensions 2, 8, and 24. 

Conjecture 7.4. The numerical method described above gives bounds 
that converge (as m — > oo) to the optimal bounds obtainable using Theorem 3.2. 



8. Uniqueness 

It is natural to ask whether the densest sphere packing in W n is unique. 
Of course, it is trivially not unique, since for example removing a single sphere 
does not change the global density. However, it is conjectured that Eg and the 
Leech lattice are unique among periodic packings. By contrast, in M 3 there are 
infinitely many distinct periodic packings of maximal density. Our techniques 
can be used to prove this uniqueness in 8 and 24 dimensions, given the following 
slight strengthening of Conjecture 7.3. At the same time, we will deal with the 
hexagonal lattice, although uniqueness in that case has long been known. 

Let A2, A§, and A24 denote the isodual scaling of the hexagonal lattice, 
the Eg root lattice, and the Leech lattice, respectively. 

Conjecture 8.1. Forn £ {2, 8, 24}, there exists a function that satisfies 
the hypotheses of Theorem 3.2 to prove that A n is the densest packing in W 1 . 
Furthermore, this function and its Fourier transform have roots only at the 
vector lengths in A n . 

First, let n be 8 or 24, so that A n is an even unimodular lattice (we will 
discuss the n = 2 case below). Let / be a function satisfying the conclusions 
of Conjecture 8.1. Suppose we have a maximally dense packing given by the 
translates of a lattice A by vectors v\, . . . , v at, whose differences are not in A. 



NEW UPPER BOUNDS ON SPHERE PACKINGS I 



707 



Without loss of generality, we can assume that |A| = N and v\ = 0. Note that 
|A| = N implies that the packing uses balls of the same radius as those in A n . 
By Poisson summation, 



l<j,k<N xeA 1 1 t&A* 



i<i<JV 



and the /(0) terms cancel the / (0) term, so we can draw the usual conclusions 
from the inequalities on each side. In particular, each vector x + v j — Vk must 
occur at a root of /. Now we can apply the following lemma: 

Lemma 8.2. Suppose S is a subset of W 1 such that G S, there are n 
linearly independent vectors in S, and for all x,y £ S, the distance \x — y\ is 
the square root of an even integer. Then the subgroup of M. n generated by S is 
an even integral lattice. 

Recall that an integral lattice is a lattice in which the inner product of 
each pair of vectors is an integer; it is even if every vector has even norm. 

Proof. For all x,y £ S, their inner product (x, y) is an integer, because 

(x,y) = (\x - 0| 2 + \y- 0| 2 - \x - y\ 2 )/2 G Z, 

and the norm \x\ 2 of any element x G S is an even integer. It follows that 
the same facts hold for integer linear combinations of elements of 5. The 
restriction on norms implies that S generates a discrete subgroup of M. n , and 
hence a lattice (because S spans R n by assumption). Now the conditions above 
on inner products and norms amount to the definition of an even integral 
lattice. □ 

Let L be the subgroup of W 1 generated by our periodic packing. By 
Lemma 8.2, it is an even integral lattice. In any integral lattice, the covolume 
is always the square root of an integer, since its square is the determinant of 
a Gram matrix, which is an integral matrix. Thus, L has at most one point 
per unit volume in W l , with equality if and only if L is unimodular. However, 
the periodic packing has one sphere per unit volume in W 1 , because |A| = N. 
It follows that the periodic packing is in fact the lattice packing determined 
by L: if any spheres from the lattice packing were missing from the periodic 
packing, then by periodicity the numbers of spheres per unit volume would be 
strictly smaller. Thus, our packing comes from an even unimodular lattice, 
and that lattice must have minimal norm 2 in M 8 and 4 in M 24 for the density 
to be right. Such lattices are unique (see Chapters 16 and 18 of [CS]). Thus, 
we have shown that Conjecture 8.1 implies that Eg and the Leech lattice are 
the only periodic packings of maximal density in M 8 and M 24 , respectively. 
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For M 2 , this argument requires a slight modification, because the isodual 
scaling of the hexagonal lattice is not an even unimodular lattice. However, the 
modification is not hard. As above, in any maximally dense periodic packing, 
the distances between points must occur among the distances in the hexagonal 
lattice. Of course, we can choose any scaling we prefer; the most convenient is 
that of the A2 root lattice, given by 

M = {(xo, xi,x 2 ) € Z 3 : x + x\ + x 2 = 0}, 

because A2 is an even integral lattice. As above, our periodic packing must 
then be contained in an even integral lattice L. However, the unimodularity 
argument no longer applies. Fortunately, we do not need it: because L is even, 
its minimal norm is at least 2, so L determines a sphere packing with spheres 
of the same radius as in our periodic packing. This sphere packing contains the 
original periodic packing. If the periodic packing did not use all these spheres, 
then its density would be lower than that of L. Thus, it is a lattice packing, 
and it is well known and easy to prove that A2 is the unique densest lattice in 
two dimensions. 

It is worth noting that the arguments above do not require the full strength 
of Conjecture 8.1. In particular, we never make use of restrictions on the 
roots of the Fourier transform. Furthermore, in the n = 2 case, we do not 
even require as strong a condition on the function's roots: A2 does not have 
vectors of every even norm, but our argument allows the function to have roots 
corresponding to the missing norms. 

Appendix A. Technicalities about density 

In this appendix, we provide precise statements and references for what it 
means for a packing to have a density and whether there is a maximally dense 
packing. Much of our discussion closely follows Section I of [K]. 

Let V be any sphere packing in W 1 . We say that V has density A if for 
all p £ K", we have 

,. vol(B(p, r) n V) 

A = hm v - 1 , — -, 

r^oo vol B(p,r) 

where B(p, r) is the ball of radius r centered at p, and B(p, r) PI V consists of 
those parts of the balls in V that lie within B(p,r). It is proved in [Gr] that 
if this limit exists for one p, then it exists for all p and is equal for all p. We 
say that the packing has uniform density if the limit exists uniformly for all p. 
In that case, [Gr] shows that for every compact set R that is the closure of its 
interior and every point p, 

,. vol((rR+p) n?) 
A = hm — , 

r->oo vol rti 

where of course rR + p denotes R scaled by a factor of r and translated by p. 
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Although not every packing has a density, every packing has an upper 
density, defined by 

vol(B(p,r)nV) 
A = iim sup sup — — r . 

r^oo pgRn vol B(p,r) 

It is proved in [Gr] that the supremum of all upper densities is achieved by a 
uniformly dense packing. 

Periodic packings are the most convenient ones for our purposes. Using the 
above results, it is easy to see that they come arbitrarily close to the greatest 
possible density, as follows. Suppose A is the maximum packing density in R n , 
and let V be a uniformly dense packing of density A. Let R be the fundamental 
parallelotope of any lattice Act". We know that 

.. voUrRnV) 
A = hm -. 

r^oo vol rR 

Let e > 0. If we choose r large enough, then the total volume of the spheres 
in V that lie entirely within rR is within e vol rR of Avolri?, because only 
a negligible fraction of the spheres can intersect the sides of rR. Define a 
periodic packing V' by taking all the spheres of V that lie entirely within rR, 
and also including all translations of them by rA. Then this periodic packing 
has density at least A — e. 

Appendix B. Other convex bodies 

Our methods are not limited to studying sphere packings, but instead 
apply to packings with translates of any convex, symmetrical body. In fact, it 
is straightforward to prove the following generalization of our main theorem: 

Theorem B.l. Let C be a convex body in R n , symmetric with respect to 
the origin. Suppose f : R n — > R is an admissible function, is not identically 
zero, and satisfies the following two conditions: 

(1) f{x) < forx^C, and 

(2) f(t) > for all t. 

Then all packings with translates of C have density bounded above by 

vol(C)/(0) 
2«/(0) 



Unfortunately, when C is not a sphere, there does not seem to be a good 
analogue of the reduction to radial functions in Theorem 3.1. That makes these 
cases somewhat less convenient to deal with. There always exist functions that 



710 



HENRY COHN AND NOAM ELKIES 



prove the trivial density bound of 1: let Xc/2 ^ e * ne characteristic function of 
the scaled body C/2, and let / = Xc/2 * ^c/2- 

Appendix C. Numerical data 

Table 3 compares the best packings known, the previous best upper bound 
known, and our bound. The packings and previous bounds are from [CS] (see 
in particular Tables 1.1(a) and 1.1(b) on pages xix and xx). Table 4 lists the 
choices of forced double roots that lead to our new bounds. 
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Dimension 


Best Packing Known 


Rogers' Bound 


New Upper Bound 


i 


U.O 


U.O 


U.O 


o 
z 


U.zooOo 


U.zoobo 


U.zoobo 


Q 

o 


n 1 7£7Q 
U. 1 1 i O 


U.lo4f 


n i q^i (\ 
U.18010 


A 

4 


n lot; 
U.lzO 


n 191 07 
U.lolz / 


n iqi or 
U.lolzO 


c 



U.Uoooy 


u.uyys i 


u.uyy i o 





n n7oi 7 
U.U / zl ( 


n HQ1 1 o 

U.Uollz 


U.UoUo4 


7 


U.UOzO 


u.uoyoi 


U.Uoyoo 


Q 
O 


U.UozO 


U.UbozO 


U.UbzOl 


o 


n n/i/ii q 

U.U44iy 


U.UOUU ( 


u.uoyuu 


1 n 

1U 


u.uoyuo 


u.uoyoo 


U.U0OU4 


11 


U.UoOlO 


U.UbloO 


u.uoyoz 


1 
1Z 


U.Uo < U4 


u.uoooy 


u.uoz / y 


lo 


U.UoOlO 


U.U ( zoo 


U.Ubo ( U 


1 A 
14 


U.UoOUo 


U.Uoz <o 


U.U m 0U 


10 


U.U44iy 


u.uy/oo 


u.usyyy 


1 R 
10 


U.UozO 


nil 77/1 
U.ll / (4 


U.1U ( oo 


1 1 


U.UOzO 


fl 1 AP.OA 
U.140Z4 


n 1Q1 
U.lolOU 


io 


U.U ( OUo 


n 1 QfiOO 

U.loOzy 


U.lbOUo 


1 Q 

iy 


U.Uoooy 


n o/i Qnc 
U.z4oUo 


n 01 ono 
U.zlzUz 


on 
zU 


n A Q~\ KA 

U.lol04 


U.OZ404 


n 07c ^ ^ 
U.Z /oOO 


zl 


n 1 7£7Q 

U.l ( ( o 


(AAA OQO 
U.44zoy 


u.o ( ooy 


00 

zz 


U.ooz04 


n R1 700 
U.01 ( zz 


n m oqi 
U.Olzol 


OQ 

zo 


U.O 


n 277K7 

U.o (lot 


n 71 RC\A 

U. ( 10U1 


OA 

z4 


1 n 
l.U 


1 070/11 

l.z ( z41 


i m one 

l.uiyyo 


ZO 


fl 7H71 1 

U. /U / 11 


1 Q70Q 

i.o <yo 


i lonm 
1.45UU1 


26 


0.57735 


2.8268 


2.18614 


27 


0.70711 


4.3252 


3.28537 


28 


1.0 


6.7295 


5.02059 


29 


0.70711 


10.642 


7.79782 


30 


1.0 


17.094 


12.30390 


31 


1.2095 


27.880 


19.71397 


32 


2.5658 


46.147 


32.06222 


33 


2.2220 


77.487 


52.90924 


34 


2.2220 


131.94 


88.55925 


35 


2.8284 


227.71 


150.29783 


36 


4.4394 


398.25 


258.54994 



Table 3. Lower and upper bounds on center density, using m = 6. 
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Dimension 


2vrr 2 


2vTZ 2 


2vrz| 


2vrz| 


2vrz| 


2irzl 


2vrz 6 2 


2 


7.25520 


21.77 


29.02 


50.79 


65.34 


90.19 




3 


8.19385 


22.00 


31.63 


49.10 


61.80 


80.01 




4 


9.10543 


22.46 


33.49 


48.69 


62.03 


79.73 




5 


9.99512 


23.04 


34.87 


48.90 


63.04 


80.44 




6 


10.86682 


23.70 


35.98 


49.43 


64.07 


81.61 




7 


11.72351 


24.41 


36.94 


50.13 


65.04 


82.94 




8 


12.56674 


25.14 


37.74 


50.50 


64.08 


79.03 


99.37 


9 


13.39945 


25.89 


38.57 


51.30 


64.99 


80.13 


100.44 


10 


14.22261 


26.66 


39.38 


52.13 


65.89 


81.24 


101.56 


11 


15.03741 


27.44 


40.18 


52.97 


66.79 


82.34 


102.71 


12 


15.84481 


28.23 


40.98 


53.82 


67.69 


83.44 


103.89 


13 


16.64565 


29.02 


41.77 


54.67 


68.61 


84.51 


105.08 


14 


17.44063 


29.82 


42.57 


55.53 


69.53 


85.58 


106.28 


15 


18.23036 


30.61 


43.37 


56.39 


70.46 


86.65 


107.48 


16 


19.01538 


31.42 


44.18 


57.25 


71.39 


87.71 


108.68 


17 


19.79617 


32.22 


44.99 


58.12 


72.34 


88.77 


109.88 


18 


20.57315 


33.02 


45.80 


58.99 


73.29 


89.84 


111.08 


19 


21.34670 


33.82 


46.62 


59.86 


74.24 


90.90 


112.28 


20 


22.11717 


34.63 


47.45 


60.74 


75.20 


91.97 


113.48 


21 


22.88488 


35.44 


48.27 


61.63 


76.17 


93.04 


114.68 


22 


23.65012 


36.24 


49.10 


62.51 


77.14 


94.12 


115.88 


23 


24.41314 


37.05 


49.94 


63.41 


78.12 


95.20 


117.08 


24 


25.17419 


37.86 


50.78 


64.30 


79.09 


96.28 


118.29 


25 


25.93349 


38.66 


51.62 


65.20 


80.08 


97.37 


119.50 


26 


26.69125 


39.47 


52.46 


66.11 


81.07 


98.46 


120.71 


27 


27.44766 


40.28 


53.31 


67.02 


82.06 


99.55 


121.92 


28 


28.20290 


41.09 


54.17 


67.93 


83.06 


100.64 


123.13 


29 


28.95711 


41.90 


55.02 


68.85 


84.06 


101.74 


124.34 


30 


29.71046 


42.72 


55.88 


69.77 


85.06 


102.84 


125.56 


31 


30.46307 


43.53 


56.74 


70.69 


86.07 


103.94 


126.77 


32 


31.21508 


44.34 


57.61 


71.62 


87.08 


105.05 


127.99 


33 


31.96659 


45.16 


58.47 


72.56 


88.09 


106.16 


129.21 


34 


32.71772 


45.97 


59.34 


73.49 


89.10 


107.27 


130.43 


35 


33.46857 


46.79 


60.21 


74.43 


90.12 


108.38 


131.65 


36 


34.21922 


47.60 


61.09 


75.37 


91.14 


109.49 


132.88 



Table 4. 



Values of 2irr 2 



and 2irz; 



,f for 1 < i < 6. 
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